Bichromatically driven double well: parametric perspective of the control landscape 
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We numerically construct and study the control landscape of a [u>, 2u>) bichromatically driven 
double well in the presence of strong fields. The control landscape is obtained by correlating the 
overlap intensities between the floquet states and an initial phase space coherent state with the 
parametric motion of the quasienergies i.e., intensity-level velocity correlator. "Walls" of no control, 
robust under variations of the relative phase between the fields, are seen in the control landscape 
and associated with multilevel interactions involving chaotic floquet states. 



o 
o 

(N 



Q 

U 
d 



> 

o 
o 



The driven double well system is an important model 
and an ideal testbed for several schemes that have been 
suggested for controlling coherent superposition of quan- 
tum states [1 H H H i, Q. Starting with the demon- 
stration of enhanced tunneling by Lin and BallentineQ 
and coherent destruction of tunneling 3] (CDT) by Gross- 
man et al. there has been a sustained interest in the 
driven double well due to its relevance in various ar- 
eas like quantum computing^, coherent control Q, and 
quantum dotsQ. Many of the proposed control scenario 
involve few, carefully chosen, levels which provide valu- 
able insights into the process. For example, the essence 
of CDT can be captured with a two-level floquet state 
analysis [l(|. However, such few-level schemes might be 
compromised due to several reasons [III EH EH, EH and 
one of them, of interest in this work, has to do with the 
effect of chaos in the underlying phase space [l2( . Typ- 
ically, control in the presence of strong fields requires 
one to go beyond two (or few) level analysis since flo- 
quet states delocalized over the chaotic region of the 
classical phase space are expected to play an impor- 
tant rolefl^. IT3L Il4j|. Examples include work on bichro- 
matically driven pendulum [13[ and a recent analysis of 
the adiabatic passage technique for coherent population 
transfer [3]. 

Despite the progress made towards understanding con- 
trol in the presence of chaos, a comprehensive under- 
standing of control in terms of the sensitivity to initial 
states, field strengths, and relative phases is still absent. 
The crucial object that one needs for such detailed un- 
derstanding is the "control landscape" since the structure 
of the control landscape, amongst other things, can shed 
light on the robustness of the control strategies [Hj]. To 
explain what we mean by a control landscape consider 
the Hamiltonian (% = 1) 
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H(x,p;t) = -—p 2 
2m 



Bx 4 - Dx 2 



xf(t) 

f(t) = Ai cos(uit) + A2 cos(2cj£ + cf>) 



(1) 



representing a driven double well. Our focus on bichro- 
matic control is motivated by extensive studies that have 
shown the utility of such fields in addition to the possibil- 
ity of the relative phase 4> providing an additional control 
knob[16l|. Owing to the periodicity of Eqjllthe dynam- 
ics of an initial state \z) can be analyzed |17| in terms of 
the floquet states {|x«)} and the associated quasiener- 



gies {E n }. The control landscape for fixed to and a given 
\z) is the two-parameter space (Ai, A2) = A exhibiting re- 
gions of control or lack thereof. An early example is the 
(Ai,cj) control plane for CDT which has been extended 
to construct the landscape for a bichromatically driven 
two-level system [l8[ . Recently [(|, the difference between 
the quasienergies of the floquet states connected to the 
bare tunneling states of Eq. [T] has been used to map the 
control landscape. 




FIG. 1: (Color online) (A) Decay time T versus the field 
strength Ai for the coherent state, \z) = \xo,po) = | — 1.5, 0), 
localized in the left regular island (indicated by an arrow in 
the inset showing the stroboscopic phase space at Ai = 6.4) 
and A2 = 0. Note the strong fluctuations in T with the 
plateau regions around Ai = 7.4, and 9.2. (B) and (C) top 
row shows the local phase space structure and survival prob- 
ability for representative points Ai = 6.4 and 9.2 respectively. 
The bottom row shows the effect of the weak control field (rel- 
ative phase 4> — 0) with strength A2 = 0.2 (B) and A2 = 0.4 
(C) on the decay of \z). 

The examples mentioned above clearly underscore the 
crucial role of avoided (exact) crossings between the flo- 
quet states. For strong driving fields the sheer number 
of such crossings, involving multiple floquet states, calls 
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for an appropriate measure to construct the control land- 
scape. Here we propose such a measure based on corre- 
lating the parametric motion of the floquet quasienergy 
levels (level velocities) with the overlap intensities p zn = 
(z|x„)| 2 . The overlap intensity-level velocity correla- 
tor was introduced by Tomsovkpjl, followed by detailed 
studies 20], as a sensitive measure of localization in phase 
space. We show, using Eq.[T]as example, that the correla- 
tor is an ideal measure for constructing and interpreting, 
perhaps quantitatively, the control landscape in the pres- 
ence of strong fields. 

In order to emphasize the role of chaos, consider the 
control problem of Eq.[T]with fixed parameter values^, @] 
m = 1, D = 10, B — 0.5 and, oj = 6.07. The double well 
supports about eight tunneling doublets and uj nearly 
corresponds to the energy separation between the ground 
and first excited states, neglecting the tunnel splittings. 
Since we are focusing on strong-field control, we choose 
Ai G [6.0, 11.0] for which the phase space shows two 
symmetry-related regular islands embedded in a chaotic 
sea. Earlier, Farrelly and Milligan have shownQ that 
control and suppression of tunneling can be achieved with 
A2 < Ai = 10. In Fig. QJA) the decay time T, in units 
of 2tt/u!, of a coherent state \z) localized in the left reg- 
ular island (cf. Fig. [IJA) inset) is shown as a function of 
Ai for A2 = 0. The decay time is infered from the first 
vanishing of the survival probability 



S(t) = |(z|-0(t))| 2 = ^2p zn p zm e~ 



■i(E„ 
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with the overlap intensities being p zn = |(z|x„}| 2 . De- 
spite the similar nature of the phase space, strong fluctu- 
ations in T reflect the complicated dynamics of \z) as de- 
termined by the participation of floquet states and their 
phase space nature. 

A natural question in the context of bichromatic con- 
trol is as follows. For a given T(Ai,A2 = 0), as in 
Fig. QJA), what choice of the 2w-field strength A2 and 
phase 6 will suppress the decay? Motivated by recent 
studies [2l[ one suspects that the local phase phase struc- 
ture about \z) is an important factor. To this end, in 
Fig. QJB) an d (C) we show the local phase space struc- 
ture near the left regular island and the associated 5(r) 
for two representative values of Ai. In case (B) a promi- 
nent 6:1 nonlinear resonance is observed. On addition of 
a weak second field A2 = 0.2 the nonlinear resonance dis- 
appears and S(t) does not decay to zero. On the other 
hand in case (C) a 9:2 resonance is observed for A2 = 
which vanishes at A2 = 0.4. However S(t) is hardly ef- 
fected although in both cases the discrete symmetry is 
broken (<^> = 0) upon addition of the 2w-field. Similar 
insensitivity to the local phase space structure occurs in 
the vicinity of regions exhibiting plateaus in Fig. HJA). 
In other words, attempts to control the decay of \z) in the 
plateau regions is bound to be difficult. Such plateau re- 
gions, near regular-chaotic avoided crossings, have been 
observed [13] in the context of chaos-assisted tunneling. 



Thus, chaotic floquet states are expected to be key for 
understanding the features in Fig.QJA). 
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FIG. 2: (Color online) (A) Survival probability S(t) of 
2} = | — 1.5,0} with Ai = 10, (j> = and varying A2 = 
(black), 0.7 (red), and 2.5 (blue). (B) The overlap intensity- 
level velocity correlator (cf. Eq. [3} is shown as a function of 
A2. C Z (A) exhibits a peak around A2 = 0.7 indicating sub- 
stantial localization and agreeing with the behaviour of S(t) 
in (A). (Bottom panel) Overlap intensities for A2 = 0,0.7 
and, 2.5 (left to right) indicate breaking of the symmetry for 
A 2 + 0. 

Several questions arise at this juncture. What is the 
precise role of the chaotic states for control? Are the 
plateaus robust for varying A2 and (jp. In order to ad- 
dress these questions we motivate the intensity-level ve- 
locity correlator [l9| and show that a simple criterion on 
the correlator highlights the regimes of control. The flo- 
quet states and quasienergies of Eq. Q] are parametrically 
dependent on the field strengths i.e., E n = E n (X) and 
IXn) = \Xn(^))- Thus the response of a floquet state 
to the 2w-field is measured by the parametric derivative, 
also known as level velocity, dE n /d\2- At the same time 
the decay dynamics of a coherent state \z) is dominated 
by floquet states that have appreciable overlap with \z) 
(cf. Eq. [2]). Therefore, floquet states \x n ) with substan- 
tial overlap p zn and large dE n jd\i are expected to be 
important in controlling the dynamics of \z). The quali- 
tative argument can be made quantitative by introducing 
the overlap intensity- level velocity correlator [n| 
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where <7 2 and a\ are the local variances of p zn and 
dE n jd\i respectively. The average in Eq. [3] is over all 
the floquet states but the dominant contributions are ex- 
pected to come from a finite number of them according 
to the qualitative argument above. In order to keep the 
level velocities to be zero centered, the mean is subtracted 
and C Z (X) is rescaled to a unitless quantity with unit 



variance, making it a true correlation coefficient [19( . In 
effect, C 2 (A) provides information on groups of states ex- 
hibiting common localization characteristics and is thus 
a quantitative measure for deviations from ergodicity. 

Utility of the correlator can be illustrated with the 
example discussed previously 0] with A x = 10 (fixed), 
A2 = 2.5 and, = 0. For the chosen parameters the 
phase space symmetry is broken and the suppression of 
tunneling is evident from Fig. [HA). However, Fig. [2fB) 
shows that C Z (A) as a function of A2 rises rapidly for 
small A2, peaks at A2 « 0.7, and varies slowly for larger 
values. Interestingly, similar values of C Z (A) at A2 = 0.7 
and A2 = 2.5 suggests comparable extent of localization 
in phase space. Indeed Fig. [DJA) shows that the tunnel- 
ing is significantly suppressed already around A 2 = 0.7. 
Note that for A2 = 2.5 the symmetry between the two 
regular islands is broken markedly resulting in the con- 
siderable suppression of tunneling seen in Fig. EJA). In 
contrast, for A2 = 0.7 the symmetry breaking is hardly 
perceptible in the phase space and yet there is a str ong 
effect on the quantum dynamics. This has been noted [131] 
before and Fig. shows that the correlator is capable of 
capturing such effects. 

We now construct the control landscape for the en- 
tire range of Ai shown in Fig.QjA). Since 2w-field is the 
control, the intensity- velocity correlator in Eq. [3] is eval- 
uated at a given point (Ai, A2) using the response of the 
floquet states with changing field strength A2. The re- 
sulting control landscape is shown in Fig.[3Jtogether with 
the decay time data of Fig.QJA) for relative phase = 0. 
The decay time plot shows regions which cannot be con- 
trolled by the 2w-field which, interestingly, correspond 
precisely to the plateaus in Fig. [TJA) . A noteworthy fea- 
ture of the control landscape is the existence of "walls" 
of no control which are characterized by C Z (X) < 0. In 
particular, such a wall at Ai 1=3 9.2 seems to be quite 
robust and persists for fairly large values of A2. Other 
such regions around Ai « 7.0, and 8.2 seem to break up 
for higher values of A2. Our analysis reveals that in the 
wall regions the dynamics of | z) is dominated by multiple 
floquet states delocalized in the phase space. An exam- 
ple of the multiple floquet contributions, as evident from 
the overlaps, and the chaotic nature of the associated 
Husimi representations Q is shown in Fig. [3J Note that 
the participation of the chaotic floquet states at A2 = 
(no control field) persists even for large A2. Hence, using 
the symmetry breaking property of the 2w-field for con- 
trol purposes is not very effective when chaotic floquet 
states are participating in the dynamics. On the other 
hand regions with C z {\) > indicate varying extent of 
control of the decay. It is interesting to note that the 
ealier analysis [BJ happens to be in a region away from 
such walls of no control! 

The effect of varying the relative phase </> on the con- 
trol landscape is shown in Fig. [4j The case 4> = ir/2 is 
shown in Fig.0jC) and one observes that C Z (A) < over 
most of the landscape which indicates very little control. 
This is consistent with the fact that the combined driving 
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FIG. 3: (Color online) The left column shows the control land- 
scape (bottom panel) compared to the decay time plot (top 
panel, solid line as in Fig. [T]A)). I n the decay time plot only 
points that are least affected upon adding the 2oj-field with 
strengths A2 = 0.2 (red circles), 0.7 (blue squares), and 2.1 
(green triangles) are shown. Note that the regions where there 
is a lack of control correspond to the plateaus with C z (A) < 0. 
The middle and right columns show the overlap intensities for 
A2 = and 2.1 respectively with Ai = 9.2 (fixed). Multiple 
state interactions involving chaotic floquet states are involved. 
The Husimi representations [7[ of the states are shown below 
the corresponding intensity plots with the axis range being 
identical to that of the inset in Fig.[ljA). See text for details. 



field is symmetric for cf> — tt/2. The control landscape for 
4> = tt/3, and 7r/4 are also shown in Fig. HJJB) and (A) 
respectively. Note that decreasing from the symmetric 
value of 7r/2 gradually leads to regions of control. How- 
ever, it is intriguing to see that the wall of no control 
seen in Fig. [3J for Ai k, 9.2 is present for all the values of 
shown here. For systems with effective h <^ 1 one ex- 
pects several such plateau regions and the effectiveness of 
the control would depend sensitively on the initial state 
and field parameters. 

In conclusion, we have constructed the control land- 
scape for the bichromatically driven double well in the 
presence of strong fields using a novel and highly sensi- 
tive measure. A simple criterion on the correlator is able 
to distinguish between regions of control and regions of 
no control. Lack of control is associated with the in- 
volvement of chaotic floquet states and the robustness of 
such regions of no control reitrates[ll| the need to excer- 
cise care in proposing finite-level control schemes in the 
presence of strong fields. We are currently exploring the 
usefulness of the correlator, applicable in the presence 
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of adiabatic fields as well, as a tool in formulating local 
control strategies. 



FIG. 4: (Color online) Control landscape (axes range as in 
Fig. O corresponding to (A) cj> = tt/4, (B) tt/3, and (C) tt/2. 
In (C) most of the landscape exhibits C Z (X) < indicating 
lack of control. Regions of control gradually emerge in (B) 
and (C) (see also Fig. [3] for (j> — 0). Note the robustness 
of the "wall" of no control around Ai ~ 9.2 (indicated by 
arrows) . 
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